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Abstract 

The production of forward jets of transverse momentum Pt{j) ~ Q and 
large momentum fraction Xj e t ^> x probes the onset of BFKL dynamics at 
HERA. A full 0{a 2 s ) calculation of the inclusive forward jet cross section is 
presented and compared to the expected BFKL cross section. The kinematical 
region populated by these events and the scale dependence of the fixed order 
perturbative QCD cross sections are discussed. 
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Deep-inelastic scattering (DIS) at HERA provides for an ideal place to probe strong 
interaction dynamics at small Bjorken x. One focus of interest has been the rise of the 
structure function Fz{x, Q 2 ) |]J for large 1/x. One would like to identify the power law 
growth, l/x ap ~ l , as predicted by the Balitsky-Fadin-Kuraev-Lipatov (BFKL) evolution 
equation. This evolution equation resums all leading a s In 1/x terms, as opposed to the more 
standard Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equation || which resums all 
leading a s In Q 2 terms. Unfortunately, the measurement of F2 in the HERA range is probably 
too inclusive to discriminate between BFKL and the conventional DGLAP dynamics [|J. 

A more sensitive test of BFKL dynamics at small x is expected from deep inelastic 
scattering with a measured forward jet (in the proton direction) and Pr{j) ~ Q 2 @- The 
idea is to study DIS events which contain an identified jet of longitudinal momentum fraction 
Xj et = Pz(jei) / E proton which is large compared to Bjorken x. 

One of the dominant Feynman graphs responsible for the parton evolution is shown in 
Fig. [I] for the particular case of a gluon initiated process. The x« denote the momentum 
fractions (relative to the incoming proton) of the incident virtual partons and pxi is the 
transverse momentum of emitted parton i. In the axial gauge, such "ladder-type" diagrams 
with strong ordering in transverse momenta, Q 2 ~ p\ n 3> ■ • • 3> PtU) 2 but only soft 
ordering for the longitudinal fraction x\ > X2 > ■ ■ ■ > x n ~ x are the source of the leading 
log Q 2 contribution. In the BFKL approximation transverse momenta are no longer ordered 
along the ladder while there is a strong ordering in the fractional momentum x n <C x n _i <C 
. . . < i! pa Xj et . When tagging a forward jet with Pr(j) ~ Q this leaves little room for 
DGLAP evolution while the condition xj et ^> x leaves BFKL evolution active. This leads 
to an enhancement of the forward jet production cross section proportional to (xj e t/x) ap ^ 1 
over the DGLAP expectation. 

Fig. [1| shows that the gluon ladders responsible for BFKL evolution first appear at 0(a^), 
i.e. for processes with at least four scattered partons in the final state. A conventional 
fixed order QCD calculation up to 0(a 2 ) contains no trace of BFKL dynamics and must 
be considered a background for its detection; one must search for an enhancement in the 
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FIG. 1. Gluon ladder diagram contributing to jet production in DIS. The position and kine- 
matics of the parton which can give rise to the forward jet is indicated. 

forward jet production cross section above the expectation for two- and three-parton final 
states. 

In this Letter we perform a full next-to-leading order (NLO) analysis of this "fixed 
order" background. Such an analysis has become possible with the implementation of QCD 
radiative corrections to dijet production in DIS in a fully flexible Monte Carlo program, 
MEPJET ||. The program integrates the cross sections for two-parton and three-parton 
final states numerically and at each phase space point the four-momenta of all final state 
particles are available. This allows to implement arbitrary experimental cuts and arbitrary 
jet definition schemes. Collinear and infrared divergences in the naive three-parton cross 
section are removed by slicing the three parton phase space into a hard region, where all 
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two-parton invariant masses squared are above a technical cutoff parameter s m i n , and the 
remaining soft and/or collinear cross section is added to the virtual contributions, yielding 
a finite two-parton cross section after splitting off initial state collinear divergences into 
universal crossing functions J7|. s m j„-independence of the resulting cross section is a 
powerful test of the integration program. For the forward jet phase space region to be 
discussed below we have verified this independence to better than 1%, which is the statistical 
accuracy of our Monte Carlo runs. 

Numerical results below will be presented both for leading order (LO) and NLO simula- 
tions. The LO 1-jet and 2-jet results employ the LO parton distributions of Gliick, Reya and 
Vogt HQ together with the one-loop formula for the strong coupling constant. At 0(a 2 s ) all 
cross sections are determined using the NLO GRV parton distribution functions f(xi,[i%) 
and the two loop formula for a s (/x^). With this procedure the 2-jet inclusive rate at NLO 
is simply given as the sum of the NLO 2-jet and the LO 3-jet exclusive cross sections. The 
value of a s is matched at the thresholds fiR = m q and the number of flavors is fixed to 
rif = 5 throughout, i.e. gluons are allowed to split into five flavors of massless quarks. 

Unless otherwise stated, both the renormalization and the factorization scales are tied 
to the sum of parton k^s in the Breit frame, 

m = Vf = 7;J2 k r(i) , (!) 
z i 

where {k^{i)) 2 = 2Ef(l — cos9 ip ). Here 9 ip is the angle between the parton and proton 
directions in the Breit frame. X^t(0 constitutes a natural scale for jet production in 
DIS [[Hj] because it interpolates between Q, in the naive parton model limit, and the sum of 
jet transverse momenta, when Q becomes negligible. 

We are interested in events with a forward jet with Pr{j) ~ Q an d Xj et ^> x and 
impose kinematical cuts which closely model the HI selection [Tl[ of such events. Jets are 



defined in the cone scheme (in the laboratory frame) with AR = 1 and \rj\ < 3.5. Here 
r] = — lntan(#/2) denotes the pseudo-rapidity of a jet. Unless noted otherwise, all jets must 
have transverse momenta of at least 4 GeV in both the laboratory and the Breit frames. 



Events are selected which contain a forward jet (denoted "j") in the pseudo-rapidity range 
1.735 < r)(j) < 2.9 (corresponding to 6.3° < 8(j) < 20°) and with transverse momentum 
P l T b U) > 5 GeV. This jet must satisfy 

Vjet = P z (j)/E p > 0.05 , (2) 
0-5<pW<4, (3) 

in the laboratory frame. The condition Xj e t ^> x is satisfied by requiring 

x < 0.004 . (4) 

Additional selection cuts are Q 2 > 8 GeV 2 , 0.1 < y < 1, an energy cut of E{1') > 11 GeV 
on the scattered lepton, and a cut on its pseudo-rapidity of —2.868 < rj{l') < —1.735 
(corresponding to 160° < 6(1') < 173.5°). The energies of the incoming electron and proton 
are set to 27.5 GeV and 820 GeV, respectively. 

Numerical results for the multi-jet cross sections are shown in Table Q. Without the 
requirement of a forward jet, the cross sections show the typical decrease with increasing jet 
multiplicity which is expected in a well-behaved QCD calculation. The 3-jet cross section 
in the last column constitutes only about 10% of the 2-jet cross section and both rates 
are sizable. The requirement of a forward jet with large longitudinal momentum fraction 
(xj et > 0.05) and restricted transverse momentum (0.5 < Pj<(j)/Q 2 < 4) severely restricts 
the available phase space. In particular one finds that the 1-jet cross section vanishes at LO, 
due to the contradicting x < 0.004 and Xj et > 0.05 requirements: this forward jet kinematics 
is impossible for one single massless parton in the final state. 

Suppose now that we had performed a full 0(a 2 ) calculation of the DIS cross section, 
which would contain 3-parton final states at tree level, 1-loop corrections to 2-parton final 
states and 2-loop corrections to 1-parton final states. These 2-loop contributions would 
vanish identically, once x <C Xj e t is imposed. The remaining 2-parton and 3-parton dif- 
ferential cross sections, however, and the cancelation of divergences between them, would 
be the same as those entering a calculation of 2-jet inclusive rates. These elements are 
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TABLE I. Cross sections for n-jet events in DIS at HERA at order a®, a s , and a 2 ,. The jet 
multiplicity includes the forward jet which, when required, must satisfy pt{]) > 5 GeV and the 
cuts of Eq. (HJ!). The transverse momenta of additional (non- forward) jets must only exceed cuts 
of 4 GeV (first and third column). This requirement is replaced by the condition > 4 GeV in 
the second column. No pj< cut is imposed in the 1-jet case at 0{oP s ) and the factorization scale is 
fixed to Q. See text for further details. 





with forward jet 
VtiV X t > 4 GeV 


k% > 4 GeV 


without forward jet 
p%,p l f b > 4 GeV 




1 jet 


pb 


pb 


8630 pb 


0(a s ): 


2 jet 


18.9 pb 


22.4 pb 


2120 pb 




1 jet inclusive 


100 pb 


100 pb 






2 jet inclusive 


83.8 pb 


98.3 pb 


2400 pb 




2 jet exclusive 


69.0 pb 


66.8 pb 


2190 pb 




3 jet 


14.8 pb 


31.5 pb 


210 pb 



already implemented in the MEPJET program which, therefore, can be used to determine 
the inclusive forward jet cross section within the cuts of Eqs. (|||-f|). At 0(a 2 s ) this cross 
section is obtained from the cross section for 2-jet inclusive events by integrating over the 
full phase space of the additional jets, without any cuts on their transverse momenta or 
pseudo-rapidities. Numerical results are shown in the third row of Table |. 

The table exhibits some other remarkable features of forward jet events: the NLO 2-jet 
inclusive cross section exceeds the LO 2-jet cross section by more than a factor four and 
the 3-jet rate at 0(a 2 s ) is about as large as the 2-jet rate at 0(a s ). These characteristics 
can be understood in terms of the kinematics of forward jet events. Kinematics puts severe 
constraints on the "recoil system" , the part of the final state in the 7-parton collision of 
Fig. which is left after taking out the forward jet. For x <C Xj e t, a high invariant mass 
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hadronic system must be produced by the photon-parton collision. For small scattering 
angle but large energy, E(j), of the forward jet this condition translates into 



Here my = y M 2 + p\ and y are the transverse mass and rapidity of the partonic recoil 
system. Eq. (|5|) implies that mx must be large, the more so the larger the ratio Xj et /x. On 
the other hand, the transverse momentum, px, of the recoil system is fixed by momentum 
conservation, pt = Iqt+PtC/)!, and the cross section is largest when the transverse momenta 
of both the virtual photon, qx, and of the forward jet are small. Thus, a large recoil 
transverse mass is most easily achieved by two or more partons which create a subsystem of 
large invariant mass, M, and some of these partons will manifest themselves as fairly hard 
hadronic jets. 

This fact is demonstrated in Fig. |2| where the transverse momentum and the pseudo- 
rapidity distributions of the recoil jet with the highest p l ^ b are shown, subject only to a 
nominal requirement of p l ^ b ,p^. > 1 GeV. Almost all forward jet events contain at least one 
second jet, with p l ^ b > 4 GeV, which typically falls into the central part of the detector. 

In the usual cone scheme final state collinear singularities are regulated by the AR 
separation cut while infrared singularities and initial state collinear emission are regulated by 
the px cut. In ^*p collisions the photon virtuality, Q 2 , eliminates any collinear singularities 
for initial state emission in the electron direction and therefore a large kx is as good a 
criterion to define a cluster of hadrons as a jet as its px- The dashed line in Fig. 0(a) 
shows the kx distribution in the Breit frame of the recoil jet candidate with the largest k^. 
Basically all forward jet events in this NLO analysis possess a recoil "jet" with k^ > 4 GeV 
and would thus be classified as 2-jet inclusive events in a variant of the cone scheme where 
the px > 4 GeV condition is replaced by a k^ > 4 GeV cut. This observation makes 
intuitively clear why we are able to calculate the 1-jet inclusive forward jet cross section 
with a program which calculates the dijet inclusive cross section at NLO: there exists a jet 
definition scheme in which all forward jet events contain at least one additional hard jet. 
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FIG. 2. Characteristics of the highest transverse momentum "jet" in the recoil system, i.e. 
excluding the forward jet. Distributions shown are (a) da/dpx in the lab frame (solid line) and 
da/dkT in the Breit frame (dashed line) and (b) the jets pseudo-rapidity distribution in the lab- 
oratory frame. All distributions are calculated at order a 2 . Jet transverse momentum cuts have 
been relaxed to p l ^ b ,p^ > 1 GeV. 

How reliable is the determination of the forward jet cross section at NLO? The impor- 
tance of higher order corrections can be estimated by studying the dependence of the cross 
section on the choice of factorization and renormalization scales, /ip and Our standard 
choice is n 2 R = up = \ (*)) • m Fig- we investigate variations of the 2-jet inclusive 

cross section when changing this scale by a factor, £, (solid lines) 

/4 = /4 = ^(E^w) 2 , (6) 

and we also consider renormalization and factorization scales which are proportional to the 
photon virtuality, Q 2 , (dashed lines) 

/4 = /4 = £<? 2 - (7) 
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FIG. 3. Scale dependence of the 2-jet inclusive cross section with a forward jet satisfying 
xj > 0.05 and pxj > 5 GeV in the lab frame (see text for additional cuts). Results are shown for 
H R = fSp = £Q 2 (dashed lines) and /i R = p? F = £(0.5 Yl ^t) 2 (solid lines), at LO (lower curves) and 
at NLO (upper curves). 

Two striking features of the forward jet cross section become apparent in this comparison. 
The large effective i-T-factor, K rs 5, was already noted in Table |I[ In addition one finds 
that the scale dependence is at least as strong at NLO as at LO. Both features are closely 
related. The smallness of the LO 2-jet compared to the NLO 2-jet inclusive cross section 
means that at least three final state partons are required to access the relevant part of 
the phase space. This three-parton cross section, however, has only been calculated at 
tree level and is subject to the typical scale uncertainties of a tree level calculation. Thus, 
even though we have performed a full 0(a 2 s ) calculation of the forward jet cross section at 
HERA, including all virtual effects, our calculation effectively only gives a LO estimate of 
this cross section and large corrections may be expected from higher order effects, like the 
gluon ladders in Fig. [I]. 

The size of these corrections may be estimated by comparing to BFKL calculations or 
to existing experimental results. The HI Collaboration has published such a measurement 
which was made during the 1993 HERA run with incident electron and proton energies of 
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E e = 26.7 GeV and E p = 820 GeV ||12|| . The acceptance cuts used for this measurement 
differed somewhat from the ones described before. Because of the lower luminosity in this 
early HERA run the Xj e t cut on the forward jet was lowered to 0.025 and defined in terms of 
the jet energy as opposed to the longitudinal momentum of the jet in the proton direction, 

x jet = E{j)/E p > 0.025 , (8) 

and the pseudo-rapidity range of the forward jet was chosen slightly larger, 1.735 < < 
2.949 (corresponding to 6° < 6(j) < 20°). Scattered electrons were selected with an energy of 
E{1') > 12 GeV and in the pseudo-rapidity range —2.794 < T](l') < —1.735 (corresponding to 
160° < 6(1') < 173°). Finally the Bjorken-x and Q 2 ranges were chosen as 0.0002 < x < 0.002 
and 5 GeV 2 < Q 2 < 100 GeV 2 . Within these cuts HI has measured cross sections of 
709 ± 42 ± 166 pb for 0.0002 < x < 0.001 and 475 ± 39 ± 110 pb for 0.001 < x < 0.002. 
These two data points, normalized to bin sizes of 0.0002, are shown as diamonds with error 
bars in Fig. |4]. Also included (dashed histogram) is a recent calculation of the BFKL cross 
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As shown before, the MEP JET program allows to calculate the full 1-jet inclusive forward 
jet cross sectionPjfor x <C Xj et . The LO result is shown as the dash-dotted histogram in Fig. |] 
and the NLO result is shown as the solid histogram. The shaded area corresponds to a scale 
variation 

/4 = /4 = £^(E 2 . (9) 

from £ = 0.1 to £ = 10, and indicates a range of "reasonable" expectations for the forward 
jet cross section at 0(a 2 ). 



While the BFKL results JTH] agree well with the HI data, the fixed order perturbative 



QCD calculations clearly fall well below the measured cross section, even when accounting 



1 We have checked that also for the kinematical region considered now almost all forward jet events 
contain at least one second jet with p l j; b > 4 GeV and > 4 GeV. 
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FIG. 4. Forward jet cross section at HERA as a function of Bjorken x within the HI acceptance 
cuts (see text). The solid (dash-dotted) histogram gives the NLO (LO) MEPJET result for 
the scale choice n 2 R = fip = £(0.5 ^r) 2 with £ = 1. The shaded area shows the uncertainty of the 
NLO prediction, corresponding to a variation of £ between 0.1 and 10. The BFKL result of Bartels 



et al. [13 1 is shown as the dashed histogram. The two data points with error bars correspond to 



the HI measurement 12 . 



for variations of the factorization and renormalization scales. The measured cross section is 
a factor 4 above the NLO expectation. The shape of the NLO prediction, on the other hand, 
is perfectly compatible with the HI results, and not very different from the BFKL curve 
in Fig. |j. At LO a marked shape difference is still observed, which can be traced directly 
to the kinematical arguments given before: according to Eq. fl5|) the transverse mass of the 
recoil system must increase proportional to Xj e t/x and this requires increased transverse 
momentum of the forward jet at LO. Thus, at LO, the expected cross section falls rapidly 
at small x, an effect which is avoided when additional partons are available in the final state 
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to balance the overall transverse momentum. 

We conclude that the existing HI data show evidence for BFKL dynamics in forward jet 
events via an enhancement in the observed forward jet cross section above NLO expectations. 
The variation of the cross section with x, on the other hand, is perfectly compatible with 
either BFKL dynamics or NLO QCD. Since MEPJET provides a full NLO prediction of the 
1-jet inclusive forward jet cross section for arbitrary cuts and jet definition schemes, more 
decisive shape tests may be possible as additional data become available. 
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